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We discuss the dynamics of the universe within the framework of Massive Graviton Dark Matter 
scenario (MGCDM) in which gravitons are geometrically treated as massive particles. In this modi- 
fied gravity theory, the main effect of the gravitons is to alter the density evolution of the cold dark 
matter component in such a way that the Universe evolves to an accelerating expanding regime, as 
presently observed. Tight constraints on the main cosmological parameters of the MGCDM model 
are derived by performing a joint likelihood analysis involving the recent supernovae type la data, 
the Cosmic Microwave Background (CMB) shift parameter and the Baryonic Acoustic Oscillations 
(BAOs) as traced by the Sloan Digital Sky Survey (SDSS) red luminous galaxies. The linear evo- 
lution of small density fluctuations is also analysed in detail. It is found that the growth factor of 
the MGCDM model is slightly different ('--^ 1 — 4%) from the one provided by the conventional flat 
ACDM cosmology. The growth rate of clustering predicted by MGCDM and ACDM models are 
confronted to the observations and the corresponding best fit values of the growth index (7) are 
also determined. By using the expectations of realistic future X-ray and Sunyaev-Zeldovich cluster 
surveys we derive the dark-matter halo mass function and the corresponding redshift distribution 
of cluster-size halos for the MGCDM model. Finally, we also show that the Hubble flow differences 
between the MGCDM and the ACDM models provide a halo redshift distribution departing signifl- 
cantly from the ones predicted by other DE models. These results suggest that the MGCDM model 
can observationally be distinguished from ACDM and also from a large number of dark energy 
models recently proposed in the literature. 

PACS numbers: 98.80.-k, 95.35.-fd, 95.36.-fx 



1. INTRODUCTION 

The high-quality cosmological observational data (e.g. 
supernovae type la, CMB, galaxy clustering, etc), accu- 
mulated during the last two decades, have enabled cos- 
mologists to gain substantial confidence that modern cos- 
mology is capable of quantitatively reproducing the de- 
tails of many observed cosmic phenomena, including the 
late time accelerating stage of the Universe. Studies by 
many authors have converged to a cosmic expansion his- 
tory involving a spatially fiat geometry and a cosmic dark 
sector formed by cold dark matter and some sort of dark 
energy, endowed with large negative pressure, in order to 
explain the observed accelerating expansion of the Uni- 
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verse Q-Q. 

In spite of that, the absence of a fundamental phys- 
ical theory, regarding the mechanism inducing the cos- 
mic acceleration, have given rise to a plethora of alter- 
native cosmological scenarios. Most are based either on 
the existence of new fields in nature (dark energy) or in 
some modification of Einstein's general relativity, with 
the present accelerating stage appearing as a sort of ge- 
ometric effect. 

The simplest dark energy candidate corresponds to a 
cosmological constant, A (see Q for reviews). In the 
standard concordance cosmological (ACDM) model, the 
overall cosmic fluid contains baryons, cold dark matter 
plus a vacuum energy. This model fits accurately the cur- 
rent observational data and it therefore provides an excel- 
lent scenario to describe the observed universe. However, 
it is well known that the concordance model suffers from, 
among others [lol |. two fundamental problems: 

(i) Fine tuning problem - the fact that the observed 
value of the vacuum energy density {p\ = Ac^/SttG ~ 
10"''^ GeV^) is more than 120 orders of magnitude below 
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the natural value estimated using quantum field theory 

m- 

(ii) Coincidence problem - the fact that the matter and 
the vacuum energy densities are of the same order just 
prior to the present epoch 12]. 

Such problems have inspired many authors to pro- 
pose alternative dark energy candidates such as A{t) cos- 
mologies, quintessence, fc— essence, vector fields, phan- 
tom dark energy, tachyons, Chaplygin gas and the list 
goes on (see |13l - [28| and references therein) . Naturally, in 
order to establish the evolution of the dark energy equa- 
tion of state (EoS), a realistic form of H{a) is required 
which should be constrained through a combination of 
independent dark energy probes. 

Nevertheless, there are other possibilities to explain 
the present accelerating stage. For instance, one may 
consider that the dynamical effects attributed to dark 
energy can be mimicked by a nonstandard gravity the- 
ory. In other words, the present accelerating stage of 
the universe can also be driven only by cold dark matter 
under a modification of the nature of gravity. Such a re- 
duction of the so-called dark sector is naturally obtained 
in the so-called f{R) gravity theories [1^ (see, however. 

On the other hand, general relativity predicts that 
gravitational waves are non-dispersive and propagate 
with the same vacuum light speed. These results lead 
to the common believe that the graviton (the "boson" 
for general relativity), must be a massless particle. How- 
ever, massive gravitons are features of some alternatives 
to general relativity as the one proposed by Visser [sij . 
Such theories have motivated many experiments and ob- 
servations in order to detect a possible dispersive behav- 
ior due to a non-zero graviton mass (see [s^] and Refs. 
there in). 

More recently, it was shown that the massive graviton 
approach proposed by Visser can be used to build realis- 
tic cosmological models that can then be tested against 
the available cosmological data [s^ . One of the main ad- 
vantages of such massive graviton cosmology is the fact 
that it contains the same number of free parameters as 
the concordance ACDM model, and, therefore, it does 
not require the introduction of any extra fields in its dy- 
namics. In this way, since the astronomical community 
is planning a variety of large observational projects in- 
tended to test and constrain the standard ACDM concor- 
dance model, as well as many of the proposed alternative 
models, it is timely and important to identify and explore 
a variety of physical mechanisms (or substances) which 
could also be responsible for the late-time acceleration of 
the Universe. 

In what follows we focus our attention to a cosmolog- 
ical model within Visser's massive graviton theory. In 
particular we discuss how to differentiate the massive 
graviton model from the concordance ACDM model. Ini- 
tially, a joint statistical analysis, involving the latest ob- 
servational data (SNIa, CMB shift parameter and BAO) 
is implemented. Secondly, we attempt to discriminate 



the MGCDM and ACDM models by computing the halo 
mass function and the corresponding redshift distribution 
of the cluster-size halos. Finally, by using future X-ray 
and SZ surveys we show that the evolution of the clus- 
ter abundances is a potential discriminator between the 
MGCDM and ACDM models. We would like to stress 
here that the abundance of collapsed structures, as a 
function of mass and redshift, is a key statistical test 
for studies of the matter distribution in the universe, 
and, more importantly, it can be accessed through ob- 
servations [36|. Indeed, the mass function of galaxy clus- 
ters has been measured based on X-ray surveys |37H39t , 
via weak and strong lensin g st udies |40l - |42| . using opti- 
cal surveys, like the SDSS [43,J3, as well as, through 
Sunayev-Zeldovich (SZ) effect [45]. In the last decade 
many authors have been involved in this kind of stud- 
ies and have found that the abundance of the collapsed 
structures is affected by the presence of a dark energy 
component (46ri58j . 

The paper is planned as follows. The basic elements 
of Visser's theory are presented in section [51 where 
we also introduce the cosmological equations for a flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) geome- 
try with massive gravitons. In section |3l a joint statisti- 
cal analysis based on SNe la, CMB and BAO is used to 
constraint the massive graviton cosmological model free 
parameter. The linear growth factor of matter pertur- 
bations is discussed in section |4l while in [SJ we discuss 
and compare the corresponding theoretical predictions 
regarding the evolution of the cluster abundances. Fi- 
nally, the main conclusions are summarized in section [51 



2. MASSIVE GRAVITONS COLD DARK 
MATTER (MGCDM) COSMOLOGY: BASIC 
EQUATIONS 

In this section we briefly present the main points of 
Visser's massive gravity approach [3l|. The full action is 
given by (in what follows h = c = 1) 
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where besides the Einstein-Hilbert Lagrangian and the 
Lagrangian of the matter flelds, we have the bi-metric 
Lagrangian; 

Anass(5,5o) = ^m/V^|(5o"')'''(3-5o)M^(.%')'(2.2) 

1 2 

x(g - ga)pu - 2 [{ga^T^g - .90)^.] 

where vrig is the graviton mass and {go)^u is a general 
flat metric. 

The field equations, which are obtained by variation of 
(|2.ip . can be written as: 
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where G^'^ is the Einstein tensor, T^'^ is the energy- 
momentum tensor for perfect fluid, and the contribution 
of the massive tensor to the field equations reads: 



{9-9^U-\{9nU{g^'r^ (2.4) 



x(5 - 5o)q/3 



Note that if one takes the Umit mg -> the standard 
Einstein field equations are recovered. 

Thus, from the construction of the Visser's theory, it 
can be classified as a bimetric theory of gravitation. This 
kind of theory was first studied by N. Rosen [32]]. In the 
Rosen's concept the metric g^^^ describes the geometry of 
the spacetime in the same way as in the context of the 
general relativity theory, and the second metric {g^j^i, 
(that Rosen denoted by ^^u) refers to the flat spacetime 
and describes the inertial forces. It is worth to mention 
that Rosen has shown that a bimetric theory satisfies the 
covariance and the equivalence principles, a fact that was 
also pointed out by Visser (for more discussion see (33|). 
In this way, in order to follow the Rosen's approach we 
have constrained the background metric to respect the 
Riemann-fiat condition, that is, R^^K.{go) = in such a 
way that we have no ambiguity on the choice of {go)fiv, it 
will always be chosen to be a flat metric, depending only 
on the particular coordinates we are dealing, of course. 

Regarding the energy-momentum conservation we will 
follow the same approach of Refs. [s^ [13] ■ Since the 
Einstein tensor satisfles the Bianchi identities V ^G^^ — 
0, the energy conservation law is expressed as: 
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In the above framework, the global dynamics of a 
flat MGCDM cosmology is driven by the following equa- 
tions^: 



BttGp = 3 - 
a 



(2.6) 



?7rGp = -2-- ( - 1 --m„V(a^-l), (2.7) 



where p is the energy density, p is the pressure and a{t) 
is the scale factor. 



^ In the present article we restrict our analysis to the flat cosmolo- 
gies in order to compare our results with those of the flat ACDM 
model that is the most accepted cosmological model as shown, 
e.g., by the WMAP7 data A generalization of the model for a 
non spatially flat cosmology will appear in a forthcoming article. 



From Eq. (|2.5p we get the evolution equation for the 
energy density, namely: 



P + 3H 



(P + P) 
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= 0, 
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where H = a/ a. By integrating the above equation for a 
matter dominated universe {p — 0) one obtains: 



po 3mg 
P(a) = ^ - 
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(2.9) 



where po is the present value of the energy density. As 
expected, in the limiting case mg — ^ all the standard 
FLRW expressions are recovered. 

Now, inserting (|2.9p in the modifled Friedmann equa- 
tion (|2.6I) we obtain the normalized Hubble parameter: 



Hi 



6H^ 



with 
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(2.10) 
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where Hq is the Hubble constant, fi™ is the matter 
density parameter (for baryons and dark matter fi^ = 
Pio/PcOj where pco = SHq^/SttG is the critical density 
parameter), and fig = f^(-^)^ is the present contribu- 
tion of the massive gravitons. It should be stressed that 
the last term of the above normalized Hubble function 
(|2.10p encodes the correction to the standard FLRW ex- 
pression. 

In general, using the FLRW equations, one can express 
the effective dark energy EoS parameter in terms of the 
normalized Hubble parameter 



U'DE(a) 



-1 



2 ^ dlnE 

3 da 



1 - nrna-^E-^{a) 



(2.12) 



After some simple algebra, it is also readily seen that 
the effective ( "geometrical" in our case) dark energy EoS 
parameter is given by (see psi. l62j): 



1 dlnSH^ 

wuE(a) = -1 - - —7, ■ 

3 ama 



(2.13) 



In our case, inserting Eq. ()2.1ip into Eq. (|2.13p it is 
straightforward to obtain a simple analytical expression 
for the geometrical dark energy EoS parameter: 



WDE(a) = -1 - ^ 



2 /7- lOa^ 



(2.14) 



It thus follows that in the cosmological context, the 
modified gravity theory as proposed by Visser can be 
treated as an additional effective fluid with EoS param- 
eter defined by (|2.14p . Note also that the current Hub- 
ble function has only two free parameters [Hq and fim), 
exactly the same number of free parameters as the con- 
ventional flat ACDM model. Naturally, the value of Hq 
is not predicted by any of the models and it is set to its 
observational value of Hq = 70 A km Mpc""'^ [5l. [63j. 
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3. LIKELIHOOD ANALYSIS 

Let us now discuss the statistical treatment of the ob- 
servational data used to constrain the MGCDM model 
presented in the previous section. 

To begin with, we consider the Constitution super- 
novae la set of Hicken et al. but in order to avoid 
possible problems related to the local bulk flow, we use a 
subset of this sample containing 366 SNe la all with red- 
shifts z > 0.02. The likelihood estimator is determined 
by a XsNia statistics: 



this case, the Xcmh function reads 



366 

XSNIa 

1=1 



(3.1) 



where Ui — {1 + Zi)^^ is the scale factor of the Universe 
at the observed redshift Zi, /i is the distance modulus 
fj, = m — M = Slogdi + 25 and is the luminosity 
distance^, ^^(a, fim) — ca^^ i/'h^v) ' from the 

likelihood analysis we find that — 0.266 ± 0.016 with 
xL(0^)/do/ ~ 446.5/365. 

In addition to the SNe la data, we also consider 
the BAO scale produced in the last scattering surface 
by the competition between the pressure of the cou- 
pled baryon-photon fluid and gravity. The resulting 
acoustic waves leave (in the course of the evolution) 
an overdensity signature at certain length scales of the 
matter distribution. Evidence of this excess was re- 
cently found in the clustering properties of SDSS galax- 
ies (see [64-66] ) and it provides a suitable "standard 
ruler" for constraining dark energy models. In this work 
we use the measurement derived by Eisenstein et al. 
[63 |. In particular, we utilize the following estimator 

2/3 



lis ^ 



SDSS data to be ^ = 
Xbao function can be written as: 



0.469 ± 0.017, where z 



measured from the 



a. 



0.35 [or 

0.75]. Therefore, the corresponding 



Xcmb(^m) 



[7^(0^) - 1.726]^ 
0.0182 



(3.3) 



It should be stressed that for CMB shift parameter, the 
contribution of the radiative component, (ilfla^"*, where 
rip ~ 4.174 X lO^^h^^) needs also to be considered 
[55. Note also that the measured CMB shift parame- 
ter is somewhat model dependent but such details of the 
models were not included in our analysis. For example, 
such is the case when massive neutrinos are included. 
The robustness of the shift parameter has been tested 
and discussed in [6^. In this case the best fit value is: 
= 0.263 ±0.03. 
The derived flm values from each individual probe ap- 
pear to be quite different, although within their mutual 
2cr uncertainty range. Therefore, in order to put tighter 
constraints on the corresponding parameter space of 
any cosmological model, the above probes are combined 
through a joint likelihood analysis'^, given by the product 
of the individual likelihoods according to: £tot{^7n) = 
'CsNia X 'Cbao X 'Ccmb, which translates in the joint 
function in an addition: Xtoti^m) = XsNia+xlAO+Xcmb- 
Now, by applying our joint statistical procedure for 
both cosmologies, we obtain the following best fit pa- 
rameters: 



• MGCDM model: n,n = 0.276 ± 0.012 with 
Xtoti^m)/dof ~ 448.5/367. Such results should 
be compared to those found by Alves et al. 
[33, namely: fl^ = 0.273 ± 0.015 with a 
xjoti^m,)/dof ~ 565.06/558. This difference must 
be probably attributed to the use of the Union2 
supernovae sample [tO] by the latter authors. 



• ACDM model: rj„ ^ 0.280 ± 0.010 with 
Xtoti^ni)/dof ~ 439.5/367, which is in good agree- 
ment with recent studies [1-11]. 



XBAo(^m) ^ 



[A(a„) - 0.469]- 
0.0172 



(3.2) 



The likelihood function peaks at flm — 0.306lQ g25. 

Finally, a very interesting geometrical probe of dark 
energy is provided by the angular scale of the sound hori- 
zon at the last scattering surface. It is encoded in the 
location of the first peak of the angular (CMB) power 
spectrum [i^, Elj , and may be defined by the quantity 
TZ = J1 a'^'E(a) ' "^^^ ^^ift parameter measured 

from the WMAP 7-years data [1] is 7^ = 1.726 ±0.019 at 
zis = 1091.36 [or ais = (1 + ^/J"^ 9.154 x 10"^]. In 




^ Since only the relative distances of the SNIa are accurate and 
not their absolute local calibration, we always marginalize with 
respect to the internally derived Hubble constant. 



^ Likelihoods arc normalized to their maximum values. In the 
present analysis we always report la uncertainties on the fitted 
parameters. Note also that the total number of data points used 
here is Ntot = 368, while the associated degrees of freedom are: 
dof= 367. Note that we sample Om £ [0.1, 1] in steps of 0.001. 
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FIG. 1: Normalized Hubble parameter as a function of red- 
shift. The solid line is the prediction of the MGCDM model. 
For comparison, the dashed line corresponds to the traditional 
ACDM model. 




a 



FIG. 2: Expansion history. In the upper panel we display the 
evolution of the dark energy effective EoS parameter. In the 
lower panel we compare the deceleration parameters of the 
MGCDM (solid line) and the concordance ACDM (dashed 
line) models. In the insert we show the relative deviation 
A{q — qa) of the two deceleration parameters. 

It should be mentioned here that using the BAO results 
of Percival et al. [65?|, does not change the previously 
presented constraints. 



4. MGCDM VERSUS ACDM COSMOLOGY 
A. The cosmic expansion history 

In Figure 1 we plot the normalized MGCDM Hubble 
function (solid line) as a function of redshift, which ap- 
pears quite different both in amplitude and shape with 
respect to the corresponding ACDM model expectations 
(dashed-line) . 

In figure 2 (upper panel), we present the evolution of 
the MGCDM effective dark energy EoS parameter. One 
can divide the evolution of the cosmic expansion history 
in different phases on the basis of the varying behavior 
of the MGCDM and ACDM models. We will investigate 
such variations in terms of the deceleration parameter, 
q{a) = — (1 -I- dlnH / dlna) , which is plotted in the lower 
panel of figure 2. In the inset plot we display the rela- 
tive deviation of the deceleration parameter, A(q — q\), 



between the two cosmological models. We can divide the 
cosmic expansion history in the following phases: 

• at early enough times a<, 0.1 the deceleration pa- 
rameters of both models are positive with g ~ ^a, 
which means that the two cosmological models pro- 
vide a similar expansion rate of the universe. Note 
that by taking the limit lim wosio.) = —5/3 for the 

a— ^0 

MGCDM model, while we always have wde = ^ 1 
for the ACDM model; 

• for 0.1 < a < 0.44 the deceleration parameters are 
both positive with q > q/^, which means that the 
cosmic expansion in the MGCDM model is more 
rapidly "decelerating" than in the ACDM case; 

• between 0.44 < a < 0.52 the deceleration parame- 
ters remain positive but q < q/^; 

• for 0.52 < a < 0.57 the traditional A model re- 
mains in the decelerated regime (^a > 0) but the 
MGCDM is starting to accelerate {q < 0); 

• for 0.57 < a < 0.94 the deceleration parameters 
are both negative and since q < qA, the MGCDM 
model provides a stronger acceleration than in 
the ACDM model (the opposite situation holds at 
0.85 <a< 0.94). 

Interestingly, prior to the present epoch (a > 0.94) the 
the deceleration parameter of the MGCDM model be- 
comes positive and when a = 1 we have wjje = 0, ie., 
the universe becomes again matter dominated, implying 
that the late time acceleration of the universe was a tran- 
sient phase which has already finished. 

From the inset panel of figure 2 it becomes clear that 
the MGCDM model reaches a maximum deviation from 
the ACDM cosmology prior to a ~ 0.75 and again at 
a ~ 1. Finally, the deceleration parameters at the present 
time are qo ~ 0.50 and goA — —0.58. If we go further to 
the future we find from Eq. (|2.14p that the state param- 
eter as well as the deceleration parameter diverges for 
a = \fljh. This value sets the turning point after which 
the universe begins to contract in the MGCDM model 
(for more details see (7T|.) 

B. The growth factor and the rate of clustering 

It is well known that for small scales (smaller than the 
horizon) the dark energy component (or "geometrical" 
dark energy) is expected to be smooth and thus it is fair 
to consider perturbations only on the matter component 
of the cosmic fiuid [tI]. This assumption leads to the 
usual equation for matter perturbations 

l„i + 2i/(5,„ - 47rGeffP™5„ = 0, (4.1) 

where the effect of "geometrical" dark energy is intro- 
duced via the expression of Gcfi = Gcfi(<) [see [73j.[73|]. 
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In the context of general relativity GeS coincides with 
the Newton's gravitational constant. Now, for any type 
of dark energy an efficient parametrization of the mat- 
ter perturbations (Sm oc D) is based on the growth rate 
/(a) = d\nD/d\na [t^, which has the following func- 
tional form: 



(4.2) 



where D{a) is the linear growth factor, ilm(a) = 
ilm.a~^ / E^(a) and 7 is the so called growth index (see 
Refs. 0, [el [zl, [73, III)- Since the growth factor of 
a pure matter universe (Einstein de-Sitter) has the form 
^EdS = ci, one has to normalize the different cosmologi- 
cal models such that D ~ a at large redshifts due to the 
dominance of the non-relativistic matter component. Us- 
ing the latter condition we can easily integrate Eq. (|4.2p 
to derive the growth factor [2^ 



D{a) = ae 



(4.3) 



In the present case we are working with a modification 
of Einstein's gravity instead of an extra fluid, in such a 
way the usual Poisson equation for the gravitational po- 
tential is modified due the presence of the mass term. 
In the simple case of the non-relativistic limit we have 
a Yukawa-like potential which accomplishes corrections 
to the Newtonian potential to scales of the order of the 
Compton wavelength of the graviton, A = ^- Using 
this kind of potential, the classic limit for the graviton 
mass obtained from solar system dynamics observations 
is Trig < 7.68 x 10~^^g [ll], but one of the most stringent 
constraints is obtained by requiring the derived dynam- 
ical properties of a galactic disk to be consistent with 
observations thereby yielding nig < 10~^^g. Now, 
by considering the best fit value obtained here for fig we 
have TTig 10~^^g, which is nearly 6 orders of magni- 
tude below to the previous bound. This value gives a 
Compton wavelength of the order of the horizon. The 
Compton wavelength can be seen as the physical length 
of graviton's perturbations which implies that these per- 
turbations play some role only close to the Hubble radius 
and thus they will be negligible at sub-horizon scales. In 
other words, Eqs. (|4.1[) . ()4.3|) are both valid also in the 
MGCDM model. 



TABLE I: Data of the growth rate of clustering [79|. The 
correspondence of the columns is as follows: redshift, observed 
growth rate and references. 



fobs 



Refs. 



0.15 
0.35 
0.55 
1.40 
3.00 



0.51 ±0.11 
0.70 ±0.18 
0.75 ±0.18 
0.90 ±0.24 
1.46 ±0.29 



[81, 82J 
[83] 

m 

[85] 
[86] 



Clearly in order to quantify the evolution of the growth 
factor we need to know the growth index. Since for the 



current graviton model there is yet no theoretically pre- 
dicted value of growth index, we attempt to provide a 
relevant value by performing a standard minimiza- 
tion procedure (described previously) between the ob- 
servationally measured growth rate (based on the 2dF 
and SDSS galaxy catalogs; see Table I; [zll) and that 
expected in the MGCDM cosmological model, according 
to: 



fobsiZi) - /model (2:^,7) 



(4.4) 



where ai is the observed growth rate uncertainty. Note 
that for comparison we perform the same analysis also 
for the ACDM model. 




FIG. 3: Upper Panel: Comparison of the observed (solid cir- 
cles _79], (see Table I) and theoretical evolution of the growth 
rate of clustering f{z). The hnes correspond to the MGCDM 
(solid curve) and the ACDM (dashed curve) models. Bot- 
tom Panel: The evolution of the growth factor, with that 
corresponding to the MGCDM model (7 — 0.56) showing a 
~ 1 — 4% difference with respect to that of the ACDM model 
(7a = 0.62), especially at large redshifts (2: > 1). Errorbars 
are plotted only for the MGCDM model in order to avoid 
confusion. 

In Figure 3 (upper panel), we present the measured 
fobs{z) (filled symbols) with the estimated growth rate 
function, f{z) = il7„(2;), for the two considered cosmolog- 
ical models. Notice, that for the MGCDM cosmological 
model (solid line) we use = 0.276 and for the A case 
(dashed line) fim — 0.280, which are the values provided 
by our likelihood analysis of section 3. In the inset panel 
of figure 3 we plot the variation of Ax^ = X^(7) — Xmin(7) 
around the best 7 fit value. For the MGCDM model we 
find 7 = 0.56+J^:i^ ix^/dof ~ 0.69), while for the ACDM 
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model we obtain 7a = 0.62ta \l {x^/dof ~ 0.75), which 
is somewhat greater, but within la, of the theoretically 
predicted value of 7a — 6/11. Such a discrepancy be- 
tween the theoretical ACDM and observationally fitted 
value of 7 has also been foun d by other authors. For ex- 
ample, Di Porto & Amendola jlj obtained 7 = O.eoto gg, 
while Nesseris & Perivolaropoulos ff^, based on mass 
fluctuations inferred from independent observations at 
different redshifts, found 7 — O-dlt-o'iT- If such a sys- 
tematic difference between the measured and the the- 
oretical 7 ACDM values is due to observational uncer- 
tainties or the method used to estimate the observed 
7, then one may expect a similar systematic difference 
to affect the measured 7 value for the MGCDM model, 
pointing to a probably more realistic value for this model 
of 7 ~ 0.49. Since however this value is within the la 
observational uncertainty, we will consider the originally 
fitted MGCDM 7 value as the nominal one. 

Using the above best fit 7 values we present, in the 
lower panel of figure 3, the growth factor evolution de- 
rived by integrating Eq. (|4.3|) for the two cosmologi- 
cal models (MGCDM-sohd and ACDM-dashed). The er- 
ror bars correspond to the la uncertainty of the fitted 
7 values. Note that the growth factors are normalized 
to unity at the present time. The difference between 
the fitted growth factors lies, at redshifts z > 1, in the 
interval ~ 1 — 4%, while when using the theoretically 
predicted ACDM value of 7a — 6/11 the difference is less 
than 1.5%. For a consistent treatment of the two models 
and for the corresponding comparison of their respective 
mass functions and halo redshift distributions we will use, 
throughout the rest of the paper, the observationally de- 
rived 7 values, ie., 7a — 0.62 and 7mgcdm — 0.56. 



5. COMPARE THE CLUSTER HALO 
ABUNDANCES 

It is important to define observational criteria that will 
enable us to distinguish between the MGCDM model and 
the concordance ACDM cosmology. An obvious choice, 
that has been extensively used, is to compare the theo- 
retically predicted cluster-size halo redshift distributions 
and to use observational cluster data to distinguish the 
models. Recently, the halo abundances predicted by a 
large variety of DE models have been coinpared with 
those corresponding to the ACDM model [la, H^]. As 
a result, such analyses suggest that many DE models ex- 
plored in this study (including some of modified gravity) 
are clearly distinguishable from the ACDM cosmology. 

We use the Press and Schecther [87J (hereafter PSc) 
formalism, based on random Gaussian fields, which de- 
termines the fraction of matter that has formed bounded 
structures as a function of redshift. Mathematical de- 



given by: 



n(M, z)dM 



p diner 
M dM 



-fpscia)dM 



(5.1) 



tails of our treatment can be found in 57[ ; here we only 
present the basic ideas. The number density of halos, 
n(M, z), with masses within the range (M, M + 6M) are 



where /psc(o') = \/W^i^c/cr) exp(— 5^/2ct^), Sc is the lin- 
early extrapolated density threshold above which struc- 
tures collapse [13, while a^(M, z) is the mass variance 
of the smoothed linear density field, extrapolated to red- 
shift z at which the halos are identified. It depends on 
the power-spectrum of density perturbations in Fourier 
space, P{k), for which we use here the CDM form ac- 
cording to [89], and the values of the baryon density 
parameter, the spectral slope and Hubble constant ac- 
cording to the recent WMAP7 results 5]. Although the 
Press-Schecther formalism was shown to provide a good 
first approximation to the halo mass function provided 
by numerical simulations, it was later found to over- 
predict/under-predict the number of low/high mass ha- 
los at the present epoch [o^, [Ml ■ More recently, a large 
number of works have provided better fitting functions 
of f{a), some of them based on a phenomenological ap- 
proach. In the present treatment, we adopt the one pro- 
posed by Reed et al. [92j . 

We remind the reader that it is traditional to 
parametrize the mass variance in terms of ag, the rms 
mass fluctuations on scales of 8 h^^ Mpc at redshift 
z = 0. 

In order to compare the mass function predictions of 
the different cosmological models, it is imperative to use 
for each model the corresponding value of Sc and ag. It 
is well known that for the usual A cosmology Sc — 1.675, 
while Weinberg & Kamionkowski [46| provide an accu- 
rate fitting formula to estimate Sc for any DE model 
with a constant equation of state parameter. Since for 
the current graviton cosmological vacuum model the ef- 
fective dark energy EoS parameter at the present time 
is w ~ which implies that the Hubble parameter is 
matter dominated it is fair to use the Einstein de-Sitter 
value Sc — 1.685 [i^. Now, in order to estimate the 
correct model erg power spectrum normalization, we use 
the formulation developed in [s^l which scales the ob- 
servationally determined erg, a value to that of any cos- 
mological model. The corresponding MGCDM value is 
f8, MGCDM — 0.828 and it is based on erg, a = 0.804 (as in- 
dicated also in Table 1), derived from an average of a va- 
riety of recent measurements (see also the corresponding 
discussion in j57,]) which are based on the WMAP7 re- 
sults [^ , on a recent cluster abundances analysis [o^ , on 
weak-lensing results [9^ and on peculiar velocities based 
analyses [95j . 

Given the halo mass function from Eq. (|5.ip wc can 
now derive an observable quantity which is the redshift 
distribution of clusters, Af{z), within some determined 
mass range, say Mi < M/Zi-^Mq < M 2 = 10^^. This 
can be estimated by integrating, in mass, the expected 
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differential halo mass function, n{M,z), according to: 



Ma 



n(M, z)dM, 



(5.2) 



Ml 



where dV/dz is the comoving volume element. In order 
to derive observationally relevant cluster redshift distri- 
butions and therefore test the possibility of discriminat- 
ing between the MGCDM and the ACDM cosmological 
models, we will use the expectations of two realistic fu- 
ture cluster surveys: 

(a) the eRDSITA satellite X-ray survey, with a flux limit 
of: /lim = 3.3 X 10"^'' ergs s~^ cm^^, at the energy band 
0.5-5 keV and covering ^ 20000 deg^ of the sky, 

(b) the South Pole Telescope SZ survey, with a limiting 
flux density at i^o — 150 GHz of /lyojim = 5 mJy and a 
sky coverage of ^ 4000 deg^. 




FIG. 4: The expected cluster redshift distribution of the 
MGCDM (solid curve) and ACDM (dashed curve) models for 
the case of two future cluster surveys (upper panels), and the 
corresponding fractional difference with respect to the refer- 
ence ACDM model (lower panels). Errorbars are 2a Poisson 
uncertainties, while the dashed lines in the lower panel bracket 
the range due to the uncertainty of the observationally fitted 
value of 7. 

To realize the predictions of the first survey we use the 
relation between halo mass and bolometric X-ray lumi- 
nosity, as a function of redshift, provided in [o^, ie: 



L{M, z) ^ 3.087 X 10 



44 



ME{z) 



1015/i-iM, 







1.554 



ergs 



(5.3) 

The limiting halo mass that can be observed at redshift z 
is then found by inserting in the above equation the lim- 
iting luminosity, given by: L = 47r(i^/ii,„c;,, with the 



luminosity distance corresponding to the redshift z and Cf, 
the band correction, necessary to convert the bolometric 
luminosity of Eg. ljOl) to the 0.5-5 keV band of eROSITA. 
We estimate this correction by assuming a Raymond- 
Smith (1977) plasma model with a metallicity of OAZq, 
a typical cluster temperature of ^ 4 keV and a Galactic 
absorption column density of uh = 10^^ cm~^. 

The predictions of the second survey can be realized 
using again the relation between limiting flux and halo 
mass from IQal: 



f 1^0, lim 



2.592 X lO^mJy 



M 



1.876 



E^^^iz) 



is the angular diameter dis- 



where dA{z) = ^^/(l + z) 
tance out to redshift z. 

In figure 4 (upper panels) we present the expected red- 
shift distributions above a limiting halo mass, which is 
Ml = Miimit = max[10"/i"iMQ, M/], with Mj corre- 
sponding to the mass related to the flux-limit at the 
different redshifts, estimated by solving Ea. (j5.3p and 
Eq. (|5.4p for M. In the lower panels we present the 
fractional difference between the MGCDM and ACDM. 
The error-bars shown correspond to 2a Poisson uncer- 
tainties, which however do not include cosmic variance 
and possible observational systematic uncertainties, that 
would further increase the relevant variance. A further 
source of uncertainty that should be taken into account 
is related to the uncertainty of the observationally de- 
rived value of 7 (see section 4). The dashed lines in 
the lower panels of Fig. 3 bracket the corresponding num- 
ber count relative model differences due to the la un- 
certainty in the value of 7, with the lower curve cor- 
responding to (7,0-8) = (0.71,0.787) and the upper to 
(7,(78) = (0.42,0.876). 

The results (see also Table II) indicate that significant 
model differences should be expected to be measured 
up to z<, 1 for the case of the eROSITA X-ray survey, 
and to much higher redshifts for the case of the South 
Pole Telescope SZ survey. What is particularly interest- 
ing is the differential difference between the ACDM and 
MGCDM models, which is negative locally {z<^ 0.3), pos- 
itive at intermediate redshifts (0.4^ z<^ 1) and negative 
again for z'^ 1. This appears to be a unique signature 
of the MGCDM model, which differentiates it from the 
behaviour of a large class of DE models (see [13) and 
makes it relatively easier to be distinguished. In Table 
II, one may see a more compact presentation of our re- 
sults including the relative fractional difference between 
the MGCDM model and the ACDM model, in character- 
istic redshift bins and for both future surveys. 



6. CONCLUSIONS 

In this work, the large and small scale dynamical prop- 
erties of a flat FLRW cold dark matter cosmology, en- 
dowed with massive gravitons (MGCDM) , were discussed 
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from an analytical and a numerical viewpoints. We find 
that the MGCDM can accommodate a "dynamic phase 
transition" from an early decelerating phase (driven only 
by cold dark matter) to a late time accelerating expan- 
sion and a subsequent recent re-deceleration phase. 

Interestingly, the Hubble function of the MGCDM 
model contains only two free parameters, namely Hq and 
rim, which is the same number of free parameters as the 
ACDM model. Performing, a joint likelihood analysis us- 
ing the current observational data (SNIa, CMB shift pa- 
rameter and BAOs), we have provided tight constraints 
on the main cosmological parameter of the MGCDM 
model, i.e., flm — 0.276 ± 0.012. We then compared 
the MGCDM scenario with the conventional flat A cos- 
mology regarding the rate of clustering as well as the 
predicted halo redshift distribution. 

The main conclusions of such a comparison are: 

• At large redshifts the amplitude of the linear per- 
turbation growth factor of the MGCDM model is 
slightly different to the A solution (at a 1 — 4% 
level) , while the observationally determined growth 



index of clustering (7 ~ 0.56) is smaller than the 
corresponding fit for the A model (7A — 0.62), al- 
though within their respective la uncertainties. 

• The shape and amplitude for the redshift distribu- 
tion of cluster-size halos predicted by the MGCDM 
model is quite different from the one of a flat 
ACDM cosmology. Such a difference depends on 
redshift and has a characteristic signature that can 
discriminate the current graviton model from other 
contender DE models in the future cluster surveys. 
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Model 


(7g 


/ 


(JA/'/A/AjeROSITA 




(6M/AfA)HPT 










Z < 0.3 


0.6 <« < 0.9 


z < 0.3 


0.6 < 2 < 0.9 


1.3 < 2 < 2 


ACDM 


0.804 


0.62 


0.00 


0.00 


0.00 


0.00 


0.00 


MGCDM 


0.831 


0.56 


-0.09 


0.15±0.01 


-0.09 


0.11 


-0.03 


MGCDM 


0.875 


0.42 


0.00 


0.25±0.01 


0.00 


0.18 


-0.01 


MGCDM 


0.789 


0.71 


-0.19 


0.06±0.01 


-0.19 


0.05 


-0.05 



TABLE II: Numerical results. The I''* column indicates the cosmological model. The 2'''* and IT'^ columns lists the corresponding 

(jg and 7 values, respectively. The remaining columns present the fractional relative difFcrcnce of the abundance of halos between 
the MGCDM and the ACDM cosmology for two future cluster surveys discussed in the text. The lower two rows show results 
corresponding to the upper and lower la range of the observational 7 value uncertainty. Errorbars are 2a Poisson uncertainties 
and are shown only if they are larger than 10~^). 



